The LSC Glitch Group : Monitoring Noise 
Transients during the fifth LIGO Science Run 
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Abstract. The LIGO Scientific Collaboration (LSC) glitch group is part of the 
LIGO detector characterization eflPort. It consists of data analysts and detector 
experts who, during and after science runs, collaborate for a better understanding 
of noise transients in the detectors. Goals of the glitch group during the fifth LIGO 
science run (S5) included (1) offline assessment of the detector data quality, with 
focus on noise transients, (2) veto recommendations for astrophysical analysis 
and (3) feedback to the commissioning team on anomalies seen in gravitational 
wave and auxiliary data channels. Other activities included the study of auto- 
correlation of triggers from burst searches, stationarity of the detector noise and 
veto studies. The group identified causes for several noise transients that triggered 
false alarms in the gravitational wave searches; the times of such transients were 
identified and vetoed from the data generating the LSC astrophysical results. 



PACS numbers: 04.80Nn, 95.55.Ym 
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1. Introduction 

The "glitch group" is one of the subgroups of the Detector Characterization Committee 
within the LIGO Scientific Collaboration (LSC). In this paper we shall use the term 
"glitch" to denote any short-duration noise transient in the gravitational wave channel 
as well as transients in auxiliary channels. Glitches produced by environmental efi^ects 
or instrumental malfunctions are a source of background for transient gravitational 
wave signals, such as unmodelled bursts or compact binary coalescences. Sufficiently 
strong glitches are also responsible for loss of lock and decreased observation time. 
The glitch group was established in 2003 to characterize noise transients in LIGO. At 
times these investigations revealed causes that could be fixed and those transients were 
eliminated. The group consists of members from the analysis groups searching for short 
duration gravitational waves from coalescing binary systems, supernovae, or other 
astrophysical systems, as well as detector experts and operators from both the LIGO 
sites at Hanford and Livingston. There is substantial interaction between the glitch 
group and other detector characterization working groups such as Calibration, Data 
Quality, Dataset Reduction, Environmental Disturbances, and Hardware Injections. This 
article will focus on the activities and findings of the glitch group during the fifth LIGO 
science run ("S5") which started in November 2005 and ended in October 2007 [1]. 
In Sect. [2] we describe the tools used by the group for the diagnoses of instrument 
artifacts. In Sect.|4]we describe some of the findings that led to elimination of identified 
problems and in better data quality. Finally, we shall conclude by highlighting some 
of the post-S5 activities (Sect.|5|. 

The glitch group provided a forum for experts in data analysis and detector 
commissioning to join forces, brainstorm and assess the performance of the LIGO 
detectors during the S5 run (4th November 2005 - 1st October 2007). 

Goals of the glitch group during the fifth LIGO science run (S5) included (1) 
offline assessment of the detector data quality, with focus on noise transients, (2) veto 
recommendations for astrophysical analysis and (3) feedback to the commissioning 
team on anomalies seen in gravitational wave and auxiliary data channels. 

The glitch group activities complemented realtime investigations and onsite 
detector troubleshooting, and provided guidance to the burst and compact binary 
coalescence (CBC) analysis groups in their veto choices. 

Members of the glitch group conducted offsite shifts, each covering 3-4 days of data 
acquisition. Results from these shifts were discussed in weekly telephone conferences. 
Highlights from these shifts were also presented each week in the run coordination 
and detector characterization teleconferences. Detailed specialized investigations were 
carried out by individual glitch group members. 

2. Goals and Methods 

Glitch group members analyzed LIGO data with several near real-time algorithms, 
with latency ranging from few minutes to a day. Some of these were run using the Data 
Monitoring Tool (DMT) environment within LIGO. The DMT is a set of algorithms 
that monitor various aspects of LIGO data quality, display status information and 
record data quality statistics. The goals of these near online algorithms ranged from 
searches for gravitational wave signals from unmodelled bursts and inspirals to studies 
of detector noise. 
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Figure 1. BurstMon 6 pixel fraction as a function of time (with one entry per 
minute). BurstMon first identifies the most significant 10 percent of time-frequency 
"pixels" using a wavelet based decomposition of the data. The pixel fraction is 
defined as the fraction of these significant pixels that can be grouped into clusters 
of two or more pixels. Ideal noise typically produces isolated pixels, resulting 
in a typical pixel fraction of 0.13. Pixel fractions much greater than this value 
indicate a high rate of glitches in the data, which preferentially produce clusters 
of significant pixels. 



Block-Normal : Block-Normal is an algorithm designed to search for short- 
duration unmodelled gravitational wave bursts. It is based on a time domain analysis 
of the data and uses a Bayesian statistics figure of merit to select candidate events [4] . 
During a glitch shift, we scanned single interferometer outliers from Block-Normal 
using the event visualization tools. 

BurstMon : BurstMon is the DMT tool for monitoring the burst detection 
performance of LIGO detectors. It is closely related to the Waveburst algorithm [5] 
that is used for untriggered gravitational wave burst searches using data from second, 
third and fourth LIGO science runs. This monitor produces 3 figures of merits: a 
measure of the rate of non-stationarity called "pixel fraction" , the real-time detector 
sensitivity to gravitational wave bursts, and noise variability in various frequency 
bands. A sample plot from the "pixel fraction" figure of merit is shown in Fig. [l] 
The term "pixel" is used to denote a time-frequency bin and the pixel fraction is the 
fraction of pixels which can be grouped into clusters of two or more pixels. It indicates 
what fraction of the time-frequency volume is affected by the non-Gaussianity of the 
detector noise. For stationary Gaussian noise, its value is equal to 0.13. However for 
real LIGO data it could be as large as 1. More details on the Burstmon figures of 
merit are provided in Ref. [B]. 

InspiralMon : Online searches for inspiralling binary compact objects were done 
using matched-filter based searches for compact object mergers between 1 and 3 
Mq using second order post-Newtonian stationary phase templates [5]. All separate 
triggers within a 15 s time- window were clustered into one set. These triggers were 
not used for the actual gravitational wave search (which uses coincidence between 
detectors, a larger template bank and several signal-based vetoes), but were very 
useful for diagnostic purposes. The signal to noise ratio of the loudest trigger was 
displayed in the control room every minute. During the glitch shifts, we examined the 
loudest 20 single interferometer CBC triggers produced each day with signal-to-noise 
ratio greater than 15. 

KleineWelle : KleineWelle [5] is a single interferometer event trigger generator. 
It is based on the dyadic wavelet decomposition of a time-series. The wavelet 
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Figure 2. KleineWelle significance versus frequency (over a one-day integrated 
period). The color-scale is tlie number of events in a given significance- frequency 
bin. 
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Figure 3. Auto-correlogram of the KleineWelle triggers. The auto-correlogram 
is a histogram of the time-difference between a given trigger and all other triggers. 
A peak in this plot is an indicator of periodicity which on a time-scale of a few 
seconds is due to enhanced microseismic noise. 



transform provides time-frequency localization of signal energy represented by the 
wavelet coefficients of the decomposition. During S5, KleineWelle analyzed in 
near-realtime (and offline) the gravitational wave channel and a variety of auxiliary 
channels for the three LIGO detectors and GEO. A variety of diagnostic plots were 
produced from these triggers. Multi-dimensional classification analysis was also done 
using these triggers JS]. Plots of the the trigger rate for a given channel with 
low and high thresholds could be produced with a graphical web-based interface. 
During the glitch shifts, we explored both double-coincident (between the 2 LIGO 
Hanford interferometers) and triple-coincident KleineWelle triggers using event 
visualization tools. We also examined various other diagnostic plots such as trigger 
auto-correlations, trigger periodicities, trigger significance as a function of frequency 
etc. Anomalous features in the auto-correlation plots are usually due to enhanced 
microseismic noise. Some of these plots are shown in Figs. [2] and [3j 

NoiseFloorMon : NoiseFloorMon is a monitor to detect slow drifts in the noise 
floor |7]. It was applied to the gravitational wave channel and to various seismic 
channels. During glitch shifts we typically looked at minute trends of threshold 
crossings and cross-correlations with seismic channels. 
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Figure 4. Time-frequency scatter plot of triggers from the QQnline analysis of 
the 4 km LIGO Hanford detector. The color indicates ranges of SNR with blue 
indicating SNRs from 5 to 10, green from 10 to 20, and red greater than 20. 
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Figure 5. Band-limited seismic noise in the 0.1 - 1 Hz band at Hanford. Such 
high microseism could cause peaks in the auto-correlogram seen in Fig. Is] 



Q Online : The Q Online pipeline is an online multi-resolution time- frequency 
search for statistically significant excess signal energy. It is equivalent to a templated 
matched filter search (in the whitened signal space), whose basis functions are 
sinusoidal Gaussians of varying central time, central frequency, and the quality factor 
Q. Details on the Q-transform are provided in Ref. |T0]. The algorithm was run online 
on data from the three LIGO detectors, VIRGO and GEO. During the glitch shifts, 
we examined the trigger trends from the QOnline pipeline (See Fig. [4]) and a scan of 
the loudest event within each hour. 

Besides the above near-online analysis, we also studied online figures of merit 
(which are usually produced in realtime in the control room) such as the effective 
distance to which LIGO is sensitive to binary neutron star inspirals, as well as 
environmental factors like wind, band-limited seismic noise, etc. See Fig. [5] for a 
plot of band- limited microseismic noise usually looked at during these shifts. 

We also gained understanding of glitch mechanisms by listening to the whitened 
versions of glitch waveforms through high-quality audio systems, taking advantage 
of the fact that our search is carried out at audio frequencies. A wide range of 
input disturbances leads to glitches with no discernible differences. In fact, each 
interferometer seems to have a characteristic glitch waveform, each a variation of a few- 
cycle oscillation near 100 Hz. Further study is exploring the exceptions to this general 
rule, including longer-duration ("more musical") tones, broad-band glitches, and 
echoes. It is hoped that these studies will give a clue about the glitch mechanism(s) , 
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Figure 6. Superposed time-series and median-normalized spectrogram as it 
appears in the Event-Display of a glitcli in the cahbration channel, after the 
application of linear predictive error filter (defined in Ref. [TOj). 



still undiagnosed. 

3. Event visualization tools 

We used two event visualization tools for a better insight into the behavior of detectors 
at any particular time of interest, that provide snapshots of the LIGO auxiliary and 
environmental channels as well as the gravitational wave channel. These are similar 
to event display tools routinely used in high energy physics experiments to depict 
the tracks of particles. Times of interest included outUers from the burst and CBC 
searches, hardware signal injections, gamma-ray-burst arrival times, environmental 
injections, etc. The tools provided insight into the behavior of the detectors at a given 
time and helped identify a few "smoking gun" causes of loud glitches, data corruption 
and sources of lock-loss. 

Event -Display : The Event-Display is also a web-based event visualization 
tool which shows the time-series and frequency spectrograms of a fixed set of channels 
along with various diagnostic information on the state of the detectors at that time, 
and output from the Paramieter Estimation code. The intensity in a given time- 
frequency bin is normalized by the median. One example of this type of specialized 
spectrogram of a glitch in the calibration channel is shown in Fig. [6] 

QSccOi : QScan is used to investigate multiple detector channels around times of 
interest. QScan produces "Q spectrogram" displays, based on the same transform 
used by the QOnline analysis. For statistically significant channels, QScan produces 
thumbnails of the time-series and "Q spectrograms" in 3 different time-windows (± 
0.5 sec, ± 2 sec, and ± 8 sec.) on a webpage. The list of channels to look at can be 
defined with a configuration file. This tool has been extensively used in the control 
room by operators and science monitors to diagnose lock-losses and is also used to look 
at the various channels in the VIRGO detector. A QScein of a glitch in the voltmeter 
channel is shown in Fig. [7] 
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Figure 7. QScan of a glitch in a voltmeter channel, after data whitening with a 
linear predictive error filter (defined in Ref. |10|V 



4. Some results from the glitch group efforts 

During S5, the glitch group provided an offline forum to explore and discuss the day- 
to-day performance of the LIGO detectors, and provide commissioners with valuable 
feedback from a data analysis perspective. More importantly, these investigations 
led to the creation of several data quality flags which are used as vetoes in the 
analysis of S5 data. The identification of noise transients is particularly important 
for the development of vetoes for the burst and CBC searches. Similar work has been 
conducted by the LSC in the past [12 [13 HH [T5] . For burst searches using the third 
LIGO science run the channel containing control signals from the power recycling 
cavity was used as a veto channel p[5l . A discusson of various data quality flags and 
veto channels used for burst searches with data from the fourth LIGO science can be 
found in Ref. |T5]. 

We provide a few examples of how our work helped commissioning efforts and 
improved data quality during the fifth LIGO science run. Although a detailed 
summary of all the investigations done by the glitch group over almost two years 
is beyond the scope of this paper, we list some of the most relevant results obtained: 

—Near the start of S5, we were able to track down causes of lock loss in the 
Livingston interferometer due to "channel hopping" , when signals intended to drive 
the amplitude of an auxiliary laser was instead injected in a channel pushing the 
detector mirrors. These are similar to the example of the calibration glitch shown in 
Fig. [6] The control and data system was fixed to monitor and prevent this artifact in 
the rest of the run. 

— We found that many coincident II1-II2 glitches were also coincident with 
events in magnetometer and voltage channels (Fig. |7]) . Our investigation of several of 
these events revealed that they coincided with circuit breaker trips, shorts, and other 
faults in high-voltage transmission lines that are connected to power substations near 
Hanford. The effect of these power grid events on the interferometer was consistent 
with what we expected from coupling of the ambient magnetic field transients to the 
permanent magnets on the test masses. The times of these and other power grid 
disruptions were flagged to prevent false alarms. 

—Starting from October 2006, LI experienced periodic glitches near the beginning 



LSC Glitch Group 



8 



Hate VI lime In the hour. Sig nificangolO 



10 20 



40 5D 



PS Hate va time in tPie lioiir. 5ia nificancoSO 



































J1 _JTJ1 





10 2 30 4 5D 60 



Figure 8. Histogram of the KleineWelle LI event rate within an hour over a one 
day period when the detector experienced hourly glitches. The top and bottom 
panels shows the rate of events with KleineWelle significance greater than 10 and 
50 respectively. 



of the hour, which were recognized as due to the digital snapshots of the various 
diagnostic information about the detector which happened once every hour. This was 
found by looking at the histograms of the KleineWelle event rate as a function of 
time within an hour as shown in Fig. [8] 

— We tracked down the cause of a few outliers in the gravitational-wave channel 
to asymmetric response in the four photodiode signals used in the optical setup. 
Dedicated monitors were written following these findings to look for similar glitches 
caused by asymmetric photodiode response throughout S5. The cause of such glitches 
was believed to be due to dust along the beam path to the photo-diodes and was 
confirmed by specalized glitch injections in the post S5 period. 

— We have also done a classification of data quality flags into four categories 
with different levels of severity. These are classifled into Category 1 (which includes 
data that won't be analyzed), Category 2 (where vetoes will be applied only in post- 
processing). Category 3 (which are advisory flags used for detection confldence) and 
Category 4 (which are advisory flags used to exert caution in case of a detection 
candidate). All these flags will be used for forthcoming burst and CBC papers using 
S5 data. 

5. Conclusions and future work 

Due to the long duration of S5 run, work is still in progress to wrap up all the S5 
related glitch group efforts. The most important S5 related task still in progress 
is the creation of data quality flags and this is being done in collaboration with 
members from the Data Quality group. Another major effort is to follow-up possible 
coincident events from burst and CBC searches to assess the data quality at the time 
of the candidate and thus their statistical significance [T7]- After S5, there were a 
few externally induced glitches and environmental injections. Some work has started 
using event visualization tools to characterize these glitches. We are also providing 
guidance to the Dataset Reduction group regarding choice of channels and sampling 
rates which need to be archived for the current Astrowatch program and future LIGO 
science runs. Thanks to the systematic effort of the glitch group, many artifacts were 
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identified and the times were flagged, producing better data quality wliicli allows for 
better astrophysical results, as well as improved confidence in any candidates that may 
be identified in the future. Many of the tools developed will be used in future runs for 
automated identification of the artifacts. The group effort has been very successful, 
and will likely continue and be improved in the future. 

Acknowledgments 

LIGO was constructed by the California Institute of Technology and Massachusetts 
Institute of Technology with funding from the National Science Foundation and 
operates under cooperative agreement PIIY-0107417. We would like to thank 
the Bonneville Power Administration, particularly Mike Viles, for providing and 
interpreting records of power grid events. We are grateful to all members of the 
LIGO Scientific Collaboration for feedback and support for our work and especially to 
all the operators and members from the commissioning team who tirelessly worked for 
two long years to ensure a high duty cycle and good sensitivity during the fifth LIGO 
science run. This paper was assigned LIGO document number LIGO-P080016-01-Z. 

References 

[1] Abbott B et al (LIGO Scientific Collaboration) 2007 "LIGO: The Laser Interferometer 

Gravitational- Wave Observatory", Preprint arXiv:0711.3041 
[2] Mukherjee S (for the LIGO Scientific Collaboration) 2007 Class. Quantum Grav. 24 S701 
[3] Zweizig J LIGO Technical Document 1999 G990079-35 
[4] McNabb J et al 2004 Class. Quantum Grav. 21 S1705 
[5] Klimenko S et al 2004 Class. Quantum Grav. 21 S1685 
[6] Klimenko S LIGO Technical Document 2005 G050192-00-Z 
[7] Mukherjee S 2006 Class. Quantum Grav. 20 S925 

[8] LIGO Scientific Collaboration Technical Document 2007 T-070109-01-Z 

[9] Chatterji S, Blackburn L, Martin G and Katsavounidis E 2004 Class. Quantum Grav. 21 S1809 
[10] Chatterji S 2005 Ph.D thesis, MIT 
[11] Sylvestre J 2002 Phys. Rev. D 66 102004 

[12] Christensen N (for the LIGO Scientific Collaboration) 2005 Class. Quantum Grav. 22 S1059 
[13] Christensen N, Shawhan P and Gonzalez G (for the LIGO Scientific Collaboration) 2004 Class. 

Quantum Grav. 20 SI 747 
[14] Cadonati L and Katsavounidis E (for the LIGO Scientific Collaboration's Bursts working group) 

2003 Class. Quantum Grav. 20 S633 
[15] Di Credico A (for the LIGO Scientific Collaboration) 2005 Class. Quantum Grav. 22, S1051 
[16] Abbott B et al (LIGO Scientific Collaboration) 2007 Class. Quantum Grav. 24 5343 
[17] Gouaty R (for the LIGO Scientific Collaboration) 2008 Preprint arXiv:0805.2412, 



